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Abstract 



The current literature on sandpile models mainly deals with the abelian sandpile 
model (ASM) and its variants. We treat a less known - but equally interesting - model, 
£S) ' namely Zhang's sandpile. This model differs in two aspects from the ASM. First, addi- 

tions are not discrete, but random amounts with a uniform distribution on an interval 
[a, b]. Second, if a site topples - which happens if the amount at that site is larger than 
a threshold value E c (which is a model parameter), then it divides its entire content in 
equal amounts among its neighbors. Zhang conjectured that in the infinite volume limit, 
this model tends to behave like the ASM in the sense that the stationary measure for 
the system in large volumes tends to be peaked narrowly around a finite set. This belief 
is supported by simulations, but so far not by analytical investigations. 

We study the stationary distribution of this model in one dimension, for several values 
of a and b. When there is only one site, exact computations are possible. Our main result 
concerns the limit as the number of sites tends to infinity, in the one-dimensional case. 
We find that the stationary distribution, in the case a > E c /2, indeed tends to that of 
the ASM (up to a scaling factor), in agreement with Zhang's conjecture. For the case 
a = 0, b = 1 we provide strong evidence that the stationary expectation tends to yl/2. 



1 Introduction and main results 

With the introduction of the sandpile model by Bak, Tang and Wiesenfeld (BTW), the notion 
of self-organized criticality was introduced, and subsequently applied to several other models 
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such as forest-fire models, and the Bak-Sneppen model for evolution. In turn, these models 
serve as a paradigm for a variety of natural phenomena in which, empirically, power laws of 
avalanche characteristics and/or correlations are found, such as the Gutenberg- Richter law 
for earthquakes. See [T2] for a extended overview. 



After the work of Dhar [2], the BTW model was later renamed 'abelian sandpile model' 
(ASM), referring to the abelian group structure of addition operators. This abelianness has 
since served as the main tool of analysis for this model. 

A less known variant of the BTW-model has been introduced by Zhang [13], where 
instead of discrete sand grains, continuous height variables are used. This lattice model is 
described informally as follows. Consider a finite subset A C 7L d . Initially, every lattice site 
i € A is given an energy < Ei < E c , where E c is the so called critical threshold, and often 
chosen to be equal to 1. Then, at each discrete time step, one adds a random amount of 
energy, uniformly distributed on some interval [a, b] C [0, .E c ], at a randomly chosen lattice 
site. If the resulting energy at this site is still below the critical value then we have arrived 
at the new configuration. If not, an avalanche is started, in which all unstable sites (that is, 
sites with energy at least E c ) 'topple' in parallel, i.e., give a fraction l/2d of their energy to 
each neighbor in A. As usual in sandpile models, upon toppling of boundary sites, energy 
is lost. As in the BTW-model, the stabilization of an unstable configuration is performed 
instantaneously, i.e., one only looks at the final stable result of the random addition. 

In his original paper, Zhang observes, based on results of numerical simulation (see also 
[6]), that for large lattices, the energy variables in the stationary state tend to concentrate 
around discrete values of energy; he calls this the emergence of energy 'quasi- units'. There- 
fore, he argues that in the thermodynamic limit, the stationary dynamics should behave as 
in the discrete ASM. However, Zhang's model is not abelian (the next configuration depends 
on the order of topplings in each avalanche; see below), and thus represents a challenge from 
the analytical point of view. There is no mentioning of this fact in [EJQ2], see however [8]; 
probably they chose the usual parallel order of topplings in simulations. 

After its introduction, a model of Zhang's type (the toppling rule is the same as Zhang's, 
but the addition is a deterministic amount larger than the critical energy) has been studied 
further in the language of dynamical systems theory in [1]. The stationary distributions 
found for this model concentrate on fractal sets. Furthermore, in these studies, emergence of 
self-organized criticality is linked to the behavior of the smallest Lyapounov exponents for 
large system sizes. From the dynamical systems point of view, Zhang's model is a non-trivial 
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example of an iterated function system, or of a coupled map lattice with strong coupling. 

In this paper we rigorously study Zhang's model in dimension d = 1 with probabilistic 
techniques, investigating uniqueness and deriving certain properties of the stationary distri- 
bution. Without loss of generality, we take E c = 1 throughout the paper. In Section [2] we 
rigorously define the model for d = 1. We show that in the particular case of d = 1 and 
stabilizing after every addition, the topplings are in fact abelian, so that the model can be 
defined without specifying the order of topplings. In that section, we also include a number 
of general properties of stationary distributions. For instance, we prove that if the number 
of sites is finite, then every stationary distribution is absolutely continuous with respect to 
Lebesgue measure on (0, 1), in contrast with the fractal distributions for the model defined 
in pQ (where the additions are deterministic). 

We then study several specific cases of Zhang's model. For each case, we prove by 
coupling that the stationary distribution is unique. In Section 01 we explicitly compute the 
stationary distribution for the model on one site, with a = 0, by reducing it to the solution 
of a delay equation [3] . 

Our main result is in Section [5l for the model with a > 1/2. We show that in the infinite 
volume limit, every one-site marginal of the stationary distribution concentrates on a non- 
random value, which is the expectation of the addition distribution (Theorem 15. 5|) . This 
supports Zhang's conjecture that in the infinite volume limit, his model tends to behave 
like the abelian sandpile. Section [5] contains a number of technical results necessary for 
proving Theorem 15.51 but which are also of independent interest. For instance, we construct 
a coupling of the so-called reduction of Zhang's model to the abelian sandpile model, and we 
prove that any initial distribution converges exponentially fast to the stationary distribution. 

In Section [6J we treat the model for [a, b] = [0,1]. We present simulations that indicate 
the emergence of quasi-units also for this case. However, since in this case there is less 
correspondence with the abelian sandpile model, we cannot fully prove this. We can prove 
that the stationary distribution is unique, and we show that if every one-site marginal of the 
stationary distribution tends to the same value in the infinite volume limit, and in addition 
if there is a certain amount of asymptotic independence, then this value is y/l/2. This value 
is consistent with our own simulations. 
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2 Model definition 



We define Zhang's model in dimension one as a discrete-time Markov process with state 
space CIn '■= [0, l)^ 1 ' 2 '"''^ C [0, oo)^ 1 ' 2 '"''^ := H^v, endowed with the usual sigma-algebra. 
We write 77, £ € f^, configurations of Zhang's model and rjj for the j'th coordinate of r\. We 
interpret rjj as the amount of energy at site j. 

By P,,, we denote the probability measure on (the usual sigma-algebra on) the path 
space for the process started in rj. Likewise we use Pj, when the process is started from a 
probability measure v on fijy, that is, with initial configuration chosen according to v. The 
configuration at time t is denoted as rj(t) and its jth component as t]j(t). 

We next describe the evolution of the process. Let < a < b < 1. At time the process 
starts in some configuration 77 £ JIat. For every t = 1, 2, . . ., the configuration r](t) is obtained 
from i](t — 1) as follows. At time t, a random amount of energy Ut, uniformly distributed 
on [a, b], is added to a uniformly chosen site Xt £ {1, . . . , N}, hence P{Xt = j) = 1/N for 
all j = 1, ... , N. We assume that Ut and Xt are independent of each other and of the past 
of the process. If, after the addition, the energies at all sites are still smaller than 1, then 
the resulting configuration is in VL^ and this is the new configuration of the process. 

If however after the addition the energy of site X t is at least 1 - such a site is called 
unstable - then this site will topple, i.e., transfer half of its energy to its left neighbor and 
the other half to its right neighbor. In case of a toppling of a boundary site, this means 
that half of the energy disappears. The resulting configuration after one toppling may still 
not be in f2jv> because a toppling may give rise to other unstable sites. Toppling continues 
until all sites have energy smaller than 1 (i.e., until all sites are stable). This final result of 
the addition is the new configuration of the process in Cl^. The entire sequence of topplings 
after one addition is called an avalanche. 

We call the above model the (N, [a, 6])-model. We use the symbol T x (£) for the result 
of toppling of site x in configuration £ G Sjv- We write A u ,x(v) f° r the result of adding an 
amount u at site x of r/, and stabilizing through topplings. 

It is not a priori clear that the process described above is well defined. By this we 
mean that it is not a priori clear that every order in which we perform the various topplings 
leads to the same final configuration r/(t). In fact, unlike in the abelian sandpile, topplings 
according to Zhang's toppling rule are not abelian in general. To give an example of non- 
abelian behavior, let TV = 2 and £ = (1.2, 1.6). Then Ti(T 2 (£)) = Ti((2,0)) = (0, 1), whereas 
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r 2 (7I(0) =T 2 ((0,2.2)) = (1.1,0). 

Despite this non-abelianness of certain topplings, we will now show that in the process 
defined above, we only encounter avalanches that consist of topplings with the abelian prop- 
erty. When restricted to a certain subset of Qn, topplings are abelian, and it turns out that 
this subset is all we use. (In particular, the example that we just gave cannot occur in our 
process.) 

Proposition 2.1. The (N, [a, b]) -model is well defined. 

Proof. We will prove in two steps that all topplings actually encountered in the process are 
abelian. To this end, we first show that in Zhang's model in dimension one, we can never 
create two adjacent unstable sites by making one addition to a stable configuration and 
toppling unstable sites in any order. 

Let &n C Hjv be the set of all (possibly unstable) configurations such that between every 
pair of unstable sites there is at least one empty site, and such that the energy of any unstable 
site is smaller than 2. It is clear that by making an addition to a stable configuration, we 
arrive in Ojv- We show that, for every configuration r\ £ fijy, the resulting configuration 
after toppling of one of the unstable sites is still in Cl^. Introducing some notation, we call 
a site j of a configuration r/ 

empty if rjj = 0, 
nonempty if rjj G (0,1), 
unstable if rjj > 1. 

An unstable site i of fj can have either two empty neighbors (first case), two nonempty 
neighbors (second case) or one nonempty and one empty (third case). 

In the first case, toppling of site i cannot create a new unstable site, since \f], L < 1, but 
% itself becomes empty. Thus, if there were unstable sites to the left and to the right of i, 
after the toppling there still is an empty site between them. 

In the second and third case, the nonempty neighbor(s) of i can become unstable. Sup- 
pose the left neighbor i — 1 becomes unstable. Directly to its right, at i, an empty site is 
created. To its left, there was either no unstable site, or first an empty site and then some- 
where an unstable site. The empty site can not have been site i — 1 itself, because to have 
become unstable it must have been nonempty. For the right neighbor the same argument 
applies. Therefore, the new configuration is still in f^. 



5 



So far, we showed that in the process of stabilization after addition to a stable configu- 
ration, only configurations in Qjv are created. In the second step of the proof we show that, 
if rj E Oat and i and j are unstable sites in rj, then 

TiiTjir,)) = Tj(Ti(rj)). (1) 

To prove this, we consider all different possibilities for x. If x is not a neighbor of either i or 
j, then toppling of i or j does not change rj x , so that ([1]) is obvious. If x is equal to i or j, or 
neighbor to only one of them, then only one of the topplings changes rj x , so that again (JTJ) 
is obvious. Finally, if j; is a neighbor of both i and j, then, since r\ € f2jv, x must be empty 
before the topplings at i and j. We then have 

so that 

1 1 

%(TjW)x = + 2^ = T j(. T i(v))x 

Therefore, also in this last case ([1]) is true. 

Having established that the topplings of two unstable sites commute, it follows that the 
final stable result after an addition is independent of the order in which we topple, and hence 
A UtX (r](t)) is well-defined; see [7], Section 2.3 for a proof of this latter fact. □ 

Remark 2.2. It will be convenient to order the topplings in so-called waves [9]. Suppose 
the addition of energy at time t takes place at site k and makes this site unstable. In the 
first wave, we topple site k and then all other sites that become unstable, but we do not 
topple site k again. After this wave only site k can possibly be unstable. If site k is unstable 
after this first wave, the second wave starts with toppling site k (for the second time) and 
then all other sites that become unstable, leaving site k alone, until we reach a configuration 
in which all sites are stable. This is the state of the process at time t. It is easy to see that 
in each wave, every site can topple at most once. 

3 Preliminaries and technicalities 

In this section, we discuss a number of technical results which are needed in the sequal, 
and which are also interesting in their own right. The section is subdivided into three 
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subsections, dealing with connections to the abelian sandpile, avalanches, and nonsingularity 
of the marginals of stationary distributions, respectively. 

3.1 Comparison with the abelian sandpile model 

We start by giving some background on the abelian sandpile model in one dimension. In the 
abelian sandpile model on a finite set A C Z, the amount of energy added is a nonrandom 
quantity: each time step one grain of sand is added to a random site. When a site is unstable, 
i.e., it contains at least two grains, it topples by transferring one grain of sand to each of its 
two neighbors (at the boundary grains are lost). The abelian addition operator is as follows: 
add a particle at site x and stabilize by toppling unstable sites, in any order. We denote 
this operator by a x : {0, 1} A — > {0, 1} A . For toppling of site x in the abelian sandpile model, 
we use the symbol T x . Abelian sandpiles have some convenient properties [2]: topplings on 
different sites commute, addition operators commute, and the stationary measure on finitely 
many sites is the uniform measure on the set of so-called recurrent configurations. Recurrent 
(or allowed) configurations are characterized by the fact that they do not contain a forbidden 
subconfiguration (FSC). A FSC is defined as a restriction of rj to a subset W of A, such that 
r) x is less than the number of neighbors of x in W, for all x. In [9], a proof can be found 
that a FSC cannot be created by an addition or by a toppling. 

In the one-dimensional case on N sites, the abelian sandpile model behaves as follows. 
Sites are either empty, containing no grains, or full, containing one grain. When an empty 
site receives a grain, it becomes full, and when a full site receives a grain, it becomes unstable. 
In the latter case, the configuration changes in the following manner. Suppose the addition 
site was x. We call the distance to the first site that is empty to the left i. If there is no 
empty site to the left, then z — 1 is the distance to the boundary, j is defined similarly, but 
now to the right. After stabilization, the sites in {x — i, . . . , x + j} n {1, . . . , N} are full, 
except for a new empty site at x — i + j. Only sites in {x — x + j} PI {1, ... , N} have 
toppled. The amount of topplings of each site is equal to the minimum of its distances to 
the endsites of the avalanche. For example, boundary sites can never topple more than once 
in an avalanche. These results follow straightforwardly from working out the avalanche. 

The recurrent configurations are those with at most one empty site; a FSC in the one- 
dimensional case is a subset of A of more than lone site, with empty sites at its boundary. 

Here is an example of how a non-recurrent state on 11 sites relaxes through topplings. An 
addition was made to the 7th site; underlined sites are the sites that topple. The topplings 
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are ordered into waves (see Remark I2.2() . In the example, the second wave starts on the 5th 
configuration: 



11011121101 -» 11011202101 -> 11012020201 -» 11020121011 -» 11101121011 
-» 11101202011 -> 11102020111 -» 11110120111 -» 11110201111 
-» 11111011111. 

To compare Zhang's model to the abelian sandpile, we label the different states of a site 
j £ {1, . . . , N} in 77 E Otv as follows: 

empty(O) if r/j = 0, 

full(l) if r/,G[±,l), 

unstable(2) if r/j > 1, 

anomalous(a) if r/j € (0, |). 

Definition 3.1. T/ie reduction of a configuration n G Qn is the configuration denoted by 
7^.(77) G {0, 1, 2, a}^ 1 '" '^ corresponding to rj by 

For general < a < 6 < 1, we have the following result. 

Proposition 3.2. For any starting configuration rj G CI, there exists a random variable 
T > with P(T < 00) = 1 such that for all t > T , r/(t) contains at most one empty 
or anomalous site. Moreover, for t > T, given that r/(t) contains an empty site xt, the 
distribution of this site is uniform on {1, . . . , N}. 

To prove this proposition, we first introduce FSC's for Zhang's model. We define a FSC 
in Zhang's model in one dimension as the restriction of n to a subset W of {1, . . . , N}, in 
such a way that 2r/j is less than the number of neighbors of j in W, for all j G W. From 
here on, we will denote the number of neighbors of j in W by deg w (j). To distinguish 
between the two models, we will from now on call the above Zhang-FSC, and the definition 
given in Section [3TT1 abelian-FSC. From the definition, it follows that a Zhang-FSC in a stable 
configuration is a restriction to a subset of more than one site, with the boundary sites either 
empty or anomalous. Note that according to this definition, a stable configuration without 
Zhang-FSC's can be equivalently described as a configuration with at most one empty or 
anomalous site. 



8 



Lemma 3.3. A Zhang-FSC cannot be created by an addition in Zhang's model. 

Proof. The proof is similar to the proof of the corresponding fact for abelian-FSC's, which 
can be found for instance in [7J, Section 5. We suppose that rj(t) does not contain an FSC, 
and an addition was made at site x. If the addition caused no toppling, then it cannot create a 
Zhang-FSC, because no site decreased its energy. Suppose therefore that the addition caused 
a toppling in x. Then for each neighbor y of x 

so that n x (r\) y > 2n y + l. Also T x {rj) x = 0, and all other sites are unchanged by the toppling. 

We will now derive a contradiction. Suppose the toppling created a Zhang-FSC, on 
a subset which we call W. It is clear that this means that x should be in W, because 
it is the only site that decreased its energy by the toppling. For all j € W, we should 
have that 2T x {rj)j < deg w (j). This means that for all neighbors y of x in W, we have 
2r] y < deg w (y) — 1, and for all other j E W we have 2<r)j < deg w (j). From these inequalities 
it follows that W \ {x} was already a Zhang-FSC before the toppling, which is not possible, 
because we supposed that rj(t) contained no Zhang-FSC. 

By the same argument, further topplings cannot create a Zhang-FSC either, and the 
proof is complete. □ 

Remark 3.4. We have not defined Zhang's model in dimension d > 1, because in that case 
the resulting configuration of stabilization through topplings is not independent of the order 
of topplings. But since the proof above only discusses the result of one toppling, Lemma [3T3l 
remains valid for any choice of order of topplings. The proof is extended simply by replacing 
the factor 2 by 2d. 

Proof of Proposition [37B . If r/ already contains at most one non-full, i.e., empty or anomalous 
site, then it contains no Zhang-FSC's, and the proposition follows. Suppose therefore that 
at some time t, rj(t) contains M(t) non-full sites, with 1 < M{t) < N. We denote the 
positions of the non-full sites of r](t) by Yi(t), i = 1, . . . , M(t), and we will show that M(t) 
is nonincreasing in t, and decreases to 1 in finite time. Note that for all 1 < i < j < M(t), 
the restriction of i](t) to {Yi(t),Yi(t) + 1, . . . , Yj(t)} is a Zhang-FSC. 

At time t + 1, we have the following two possibilities. Either the addition causes no 
avalanche, in that case M(t + 1) < M(t), or it causes an avalanche. We will call the set of 
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sites that change in an avalanche (that is, all sites that topple at least once, together with 
their neighbors) the range of the avalanche. We first show that if the range at time t + 1 
contains a site y € {Yi(t), . . . , Yi + i(t)} for some i, then M(t + 1) < M(t). 

Suppose there is such a site. Then, since {Yi(t) + 1, . . . , 5^+1 (t) — 1} contains only full 
sites, all sites in this subset will topple and after stabilization of this subset, it will not 
contain a Zhang-FSC. In other words, in this subset at most one non-full site is created. But 
since Yi{t) and Y^+i(t) received energy from a toppling neighbor, they are no longer empty 
or anomalous. Therefore, M(t + 1) < M(t). 

If there is no such site, then the range is either {l,Yi(t)}, or {Y M u\ (t), N}. With the 
same reasoning as above, we can conclude that in these cases, Y\{t + 1) < resp. 
Y M(t) (t + l)>Y M(t) (t). 

Thus, M(t) strictly decreases at every time step where an avalanche contains topplings 
between two non-full sites. As long as there are at least two non-full sites, such an avalanche 
must occur eventually. We cannot make infinitely many additions without causing topplings, 
and we cannot infinitely many times cause an avalanche at x < Yi(t) or x > Y^^(t) without 
decreasing M(t), since after each such an avalanche, these non-full sites 'move' closer to the 
boundary. □ 

In the case that a > 1/2, we can further specify some characteristics of the model. We prove 
that for any initial configuration, after at most N(N — 1) time steps there is at most one 
empty site, and the other sites are is full, i.e., there are no anomalous sites. We will call 
such configurations regular. 

Proposition 3.5. Suppose a > \. Then 

1. for any initial configuration rj, for all t > N(N — 1), rj(t) is regular, 

2. for every stationary distribution fx, and for all i € {1, ... , N}, 

KVi = 0) = j^. 



In words, this proposition states that if a > 1/2, then every stationary distribution 
concentrates on regular configurations. Moreover, the stationary probability that a certain 
site i is empty, does not depend on i. Note that as a consequence, the stationary probability 
that all sites are full, is also jttt- 
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To prove this proposition, we need the following lemma. In words, it states that if 
a > 1/2 and rj contains no anomalous sites, then the reduction of Zhang's model (according 
to Definition I3.1|) behaves just as the abelian sandpile model. 

Lemma 3.6. For all u G [5, 1), for all r\ G £In which do not contain anomalous sites, and 
for allx G {1, . . . , N} 

K(AuM) = a x {n{rj)), (3) 

where a x is the addition operator of the abelian sandpile model. In both avalanches, corre- 
sponding sites topple the same number of times. 

Proof. Under the conditions of the lemma, site x can be either full or empty. If x is empty, 
then upon the addition of u > ^ it becomes full. No topplings follow, so that in that case 
we directly have H{A u , x (r])) = a x (TZ(r])). 

If 77 is such that site x is full, then upon addition it becomes unstable. We call the con- 
figuration after addition, but before any topplings 77. To check if in that case TZ(A UiX (r])) = 
a x (lZ(ri)), we only need to prove TZ(T x (fj)) = T x (TZ(f])), with 1Z(fj) x = 2, since we already 
know that in both models, the final configuration after one addition is independent of the 
order of topplings. 

In T x (fj), site x will be empty. This corresponds to the abelian toppling, because site 
x contained two grains after the addition, and by toppling it gave one to each neighbor. 
In T x (rj), the energy of the neighbors of x is their energy in rj, plus at least ^. Thus the 
neighbors of site x will in T x (fj) be full if they were empty, or unstable if they were full. Both 
correspond to the abelian toppling, where the neighbors of x received one grain. □ 

Proof of Proposition [375[ To prove part (1), we note that any amount of energy that a site can 
receive during the process, i.e., either an addition or half the content of an unstable neighbor, 
is at least 1/2. Thus, anomalous sites can not be created in the process. Anomalous sites 
can however disappear, either by receiving an addition, or, as we have seen in the proof of 
Proposition 13.21 when they are in the range of an avalanche. 

When we make an addition of at least 1/2 to a configuration with more than one non-full 
site, then either the number of non-full sites strictly decreases, or one of the outer non-full 
sites moves at least one step closer to the boundary. We note that rj contains at most N 
non-full sites, and the distance to the boundary is at most N — 1. When finally there is only 
one non-full site, then in the next time step it must either become full or be in the range of 
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an avalanche. Thus, there is a random time T' < N(N — 1) such that ry(T') is regular for 
the first time, and as anomalous sites cannot be created, by Proposition 13.21 n(t) is regular 
for all t > T'. 

For t > T' , r](t) satisfies the condition of Lemma 13.61 This means that the stationary 
distribution of the reduction of Zhang's model must coincide with that of the abelian sandpile 
model. As we mentioned in Section \'6.1\ this is the uniform measure on all configurations 
with at most one empty site. This proves part (2). □ 

3.2 Avalanches in Zhang's model 

We next describe in full detail the effect of an avalanche, started by an addition to a con- 
figuration rj(t) in Zhang's model. Let C(t + 1) be the range of this avalanche. Recall that 
we defined the range of an avalanche as the set of sites that change their energy at least 
once in the course of the avalanche (that is, all sites that topple at least once, together with 
their neighbors). We denote by T(t + 1) the collection of sites that topple at least once in 
the avalanche. Finally, C'(t + 1) C C(t + 1) denotes the collection of anomalous sites that 
change, but do not topple in the avalanche. 

During the avalanche, the energies of sites in the range, as well as Ut+\, get redistributed 
through topplings in a rather complicated manner. By decomposing the avalanche into waves 
(see Remark I2.2D . we prove the following properties of this redistribution. 

Proposition 3.7. Suppose an avalanche is started by an addition at site x to configuration 
n(t). For all sites j in C(t + 1), there exist Fij = Fij(r](t),x, Ut+i) such that we can write 



+ 1) = FijViit) + FxjUt+i + Vj(t)^jec'(t+i) 



(4) 



ier(t+i) 



with 



1. 




(5) 



ieT(t+l) 



2. for all j e C(t + 1) such that rjj(t + 1) / 0, 
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3. for all j G C(t + 1) such that r]j(t + 1) ^ 0, j > x, we have 

and similarly, F x j-\ < F x j for j < x. 

In words, we can write the new energy of each site in the range of the avalanche at time 
t + 1 as a linear combination of energies at time t and the addition f/t+i, in such a way that 
the prefactors sum up to 1 or 0. Furthermore, every site in the range receives a positive 
fraction of at least 2 T 3 ^/ 2 ! Q f the addition. These received fractions are such that larger 
fractions are found closer to the addition site. We will need this last property in the proof 
of Theorem 15.51 



Proof of Proposition {3~77\ We start with part (1). First, we decompose the avalanche started 



at site x into waves. We index the waves with ft = 1, . . . ,K, and write out explicitly the 
configuration after wave ft, in terms of the configuration after wave ft — 1. The energy of site i 
after wave k is denoted by fji t i~; we use the tilde to emphasize that these energies are not really 
encountered in the process. We define r^o = Viif) + Ut+ili= x ', note that rji t K = rji(t + 1). 

In each wave, all participating sites topple only once. We call the outermost sites that 
toppled in wave k, the endsites of this wave, and we denote them by and M' k , with 
Mk > M' k . For the first wave, this is either a boundary site, or the site next to the empty 
or anomalous site that stops the wave. Thus, M\ and M{ depend on r/(t), x and Ut+i- All 
further waves are stopped by the empty sites that were created when the endsites of the 
previous wave toppled, so that for each ft, Mk+i = — 1 and M' k+l = M' k + 1. In every 
wave but the last, site x becomes again unstable. Only in the last wave, x is an endsite, so 
that at most one of its neighbors topples. 

In wave ft, first site x topples, transferring half its energy, that is, ^fj x ^-i, to each 
neighbor. Then, if x is not an endsite, both its neighbors topple, transferring half of their 
current energy, that is, hvx±l,k-l + jVx,k-i, to their respective neighbors. Site x is then 
again unstable, but it does not topple again in this wave. Thus, the topplings propagate 
away from x in both directions, until the endsites are reached. Every toppling site in its turn 
transfers half its current energy, including the energy received from its toppling neighbor, to 
both its neighbors. Writing out all topplings leads to the following expression, for all sites 
i > x. A similar expression gives the updated energies for the sites with i < x. Note that 
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for every k > 1, rjM k +i,k-i = 0. Only when k = 1, it can be the case that site M\ + 1 was 
anomalous, so that 77^1+1,0 > 0. 



Vx,k 
VM k ,k 



1 



1 



1 



1 



i+1 j 

X] 9 i+2-n ^ fc - 1 ' f ° r i = X+l,...,M k -l, 



n=0 

0, 



(6) 



VM k -i,k + VM h +i,k-i if M k >x + 2, 

\fjx+l,k-l + \Vx,k-l + fjx+2,k-l M k = x + 1, 



I ivx,k-i + fjx+i,k-i M k = x. 

We write for all j £ C(t + 1), with fij(k) implicitly defined by the coefficients in ([6]), 

= X fij(k)fji,k-i- 
i6T(t+l) 

Since we made an addition to a stable configuration, we only encounter configurations in 
Qn- From a case by case analysis of ([6]), we claim that for all j € C(t + 1) we have 

K(Vj,k) = ^2 fijityKifjit-i); (7) 

ieT(t+l) 

the reader can verify this for all cases. 

To prove the proposition, we start with j £ C'(i + 1), for which we have 

rjj(t + 1) = fj jt i = fiji 1 )^ + Vj(t), 

ieT(t+l) 

which is dH) with i^- = fij(l). We also have, according to (J7D, 

^•(t + i)) = X fijO-Wrkfl) 

ieT(t+i) 

ier(t+i) 

since a site in + 1) becomes full in the avalanche. This proves part © of the proposition 
for such sites. 
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For all other sites in C(t + 1), we use induction in k. For wave k — 1, we make the 
induction hypothesis that 



with 



%,fc_l = ^(fc - l)*7m(*) + - l)f7t + i, 

meT((+l) 

^ F m j(k - 1) + ^(fc - 1) = 7e(% lfc -i). 

meT(m) 



(8) 



(9) 



For we then obtain 

fjj,k = 



Y fij( k )Vi,k-l 

ieT(t+i) 

Y Y F mi(k-l)f ij (k)v m (t)+ Y fij(k)F xi (k-l)U t+1 . 

meT (t+1) ieT (t+1) ieT (t+1) 



We also have 



Y fv( k ) n (Vj,k-l) 
ieT(t+l) 



Y M k ) 

ieT(t+i) 



^ ^ Fmi ( 

ieT(t+l) 



Hence, if we define 



F mj (k)= Y fij(k)F mi (k-l), 
ieT (t+1) 

then ([8]) and ([9]) are also true for wave k. For A; — 1 = 0, the hypothesis is also true, with 
F m i(0) = \ m =%- If we define Fij := Fij(K), then the first part of the proposition follows. 

To prove part (2) of the proposition, we derive a lower bound for F x j. The number K of 
waves in an avalanche is equal to the minimum of the distance to the end sites, leading to 
the upper bound K < [N/2] . 

After the first wave, © gives for all nonempty j ^ x, F x j(l) > (^) N+1 . At the start of 
the next wave, the fraction of Ut+\ present at x is equal to F xx (l) = o- Hence, after the 
second wave, even if we ignore all fractions of Ut on sites other than x, then we still have, 
again by ©, F x j{2) > \{\) N+1 • So if before each wave we always ignore all fractions of Ut 
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on sites other than x, and if we assume the maximum number of waves, then we arrive at a 
lower bound for nonempty sites j: 

F xj > (^ Nm ~Hl) N+1 > rw. 

To prove part (3) of the theorem (we only discuss the case j > x, since by symmetry, the 
case j < x is similar), we show that for every k 6 {1, ... , K}, 

F xx {k) > F x>x+1 (k) >■ ■ ■> F x>Mk -i(k) = F x>Mk +i(k), (10) 

and 

F x , Mk +l(k) > F x>Mk _ 1+ i(k - 1). (11) 

This is sufficient to prove the theorem, since for every k, site M& + 1 does not change anymore 
in the waves k + 1, . . . ,K. 

After the first wave, we have from (O that 



F XX {1) > F XsX+1 (l) > ■■■> F xMk -i(l) = F xMk+1 (l) 



so that pOj) and pij) are satisfied after the first wave. For all other waves except the last 
wave, we apply induction. 

Assume that after wave k — 1, for k < K, 

^F xx (k - 1) > F x , x+1 (k - 1) > • • • > F xMk ^-i{k - 1) = F xMk _ 1+ i(k - 1). (12) 

We have seen that this hypothesis is true after the first wave. We rewrite ©, for every 
k < K (so that > 1 and F X)X+ \{k — 1) = F XjX _i(k — 1)), as follows: 

F xx (k) = F x , x+1 (k-l) + -F xx (k-l), 
F XjX+ x(k) = -F XjX+ 2(k - 1) + -F xx (k), 

F XiX+i (k) = -F XjX+i+1 (k - 1) + - -F-r^-i-i (k) for i = 2, . . . ,M k - 1, 
F xMk (k) = 0, 
F x ,M k +i(k) = F X)Mk -i(k). 

(13) 
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From (|13p and (|12p . we find the following inequalities, each one following from the previous 
one: 

F xx (k) = F x , x+1 (k - 1) + -F xx (k - 1) < -F xx {k - 1) + -F xx (k - 1) = F xx {k - 1), 

^x,x+l(*) = ^F X}X+2 (k - 1) + -F M (fc) < -F XtX+1 (k - 1) + -^(fc - 1) = -i^fc), 

1 11 1 

F x , x+2 (k) = -F x>x+3 (k - 1) + -F XtX+1 (k) < -F X)X+2 (k - 1) + -^(fc) = F XiX+1 (k). 

For all ? = 2, . . . , M k - 1, if i^^+^/c) < i^ >x -+i-i (^) then 

i^x+i+iW = ^F x>x+i+2 {k-l) + -F x>x+ i{k) < -F XtX+ i+i(k-l) + -F XjX+ i-i(k) = F XjX+i (k). 

(14) 

Since F xx +i(k) < F x x +i—i (k) is true for i = 2, (|12j) follows for wave A;, and (jlOp is proven 
for every k < K. Moreover, we have 

1 1 

F x ,M k +l(k) = F XjMk -i(k) = -F XiMk (k - 1) + ~F Xj M k -2(.k). 

With the above derived F X} M k -i(k) < F X} M k -2(k), it follows that F Xt M k (k — l) < F Xj M k -2(k), 
so that 

F X) M k -i(k) > F XjMk {k - 1) = F x>Mk ^+i{k - 1), 

which is (jlip . 

Finally we discuss the last wave. In case Mk = we have 

F XX (K) = 0, 
F X>X+1 {K) = ^F xx (K-l) = F x , x+2 (K-l). 



In case Mk = 1 we have 

F XX (K) = If x , x+1 (K-1) + ~F xx (K-1), 
F XiX+1 (K) = 0, 

F XjX+2 (K) = F XX (K) = ^F XjX+1 (K - 1) + ^F XX {K - 1) > F XjX+1 (K - 1). 
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For all Mk > 1 we have 



F X , X+1 (K) 



F XX (K) 



^F X;X+1 (K-1) + ^F XX (K-1), 
^F X}X+2 (K - 1) + X -F XX {K - 1) < F XX (K). 



We can now use (|14p as above for i = 1, . . . , Mk — 1, so that (jlOp and (jlip follow. 



□ 



3.3 Absolute continuity of one-site marginals of stationary distributions 

Consider a one-site marginal Uj of any stationary distribution v of Zhang's sandpile model. 
It is easy to see that Vj will have an atom at 0, because after each avalanche there remains 
at least one empty site. It is intuitively clear that there can be no other atoms: by only 
making uniformly distributed additions, it seems impossible to create further atoms. Here 
we prove the stronger statement that the one-site marginals of any stationary distribution 
are absolutely continuous with respect to Lebesgue measure on (0, 1). 

Theorem 3.8. Let v be a stationary distribution for Zhang's model on N sites. Every 
one-site marginal of v is on (0, 1) absolutely continuous with respect to Lebesgue measure. 

Proof. Let A C (0, 1) be so that X(A) = 0, where A denotes Lebesgue measure. We pick a 
starting configuration rj according to v. We define a stopping time r as the first time t such 
that all non-zero energies r]i(t) contain a nonzero contribution of at least one of the added 
amounts U%, U 2 , ■ ■ ■ , Up We then write, for an arbitrary nonzero site j, 



The second term at the right hand side tends to as t — > 00 by item 2 of proposition 13.71 
We claim that the first term at the right hand side is equal to zero. To this end, we first 
observe that T)j(t) is built up of fractions of ^(0) and the additions Ui, U2, ■ ■ ■ , Up These 
fractions are random variables themselves, and we can bound this term by 



where Z± represents the (random) fraction of ryi(0) in r]j(t), and Y s represents the (random) 
fraction of U s in rjj(t). 



F u ( Vj (t) ei)< Pu(Vj(t) £A,t< t)+F v (t < t). 



(15) 




(16) 
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We clearly have that the U s are all independent of each other and of r/i(0) for all i. 
However, the U s are not necessarily independent of the Z% and the Y s , since the numerical 
value of the U s effects the relevant fractions. Also, we know from the analysis in the previous 
subsection that the Z% and Y s can only take values in a countable set. Summing over all 
elements in this set, we rewrite (1161) as 



/ n t 

Kl^2z l r H (o) + Y,y s u s eA,z i = z i ,Y s 



which is at most 

/ N 



,X N ). 



2 p„ I J2 z M°)+J2y sUs eA \ 

which, by the independence of the U s and the 77,(0), is equal to 

Klj2ziXi + J2y s U s eAjdu(xi, 

Since Y? s =i Us > 0) U s are independent uniforms, and by assumption X(A) = 0, the proba- 
bilities inside the integral are clearly zero. Since the left hand side of (|15[) is equal to Vj(A) 
for all t, we now take the limit t — > 00 on both sides, and we conclude that Vj(A) =0. □ 

Remark 3.9. The same proof shows that for every stationary measure u, and for every 
ii,-.-,ik £ {1] • • • ; -W}, conditional on i\, . . . being nonempty, the joint distribution of 
,rji k under v is absolutely continuous with respect to Lebesgue measure on (0, l) k . 



4 The (1, [a,6])-model 

In this section we consider the simplest version of Zhang's model: the (1, [a, 6])-model. In 
words: there is only one site and we add amounts of energy that are uniformly distributed 
on the interval [a, b], with < a < b < 1. 

4.1 Uniqueness of the stationary distribution 

Before turning to the particular case a = 0, we prove uniqueness of the stationary distribution 
for all [0,6] C [0,1]. 



19 



Theorem 4.1. (a) The (1, [a, b]) model has a unique stationary distribution p = p . For 
every initial distribution ¥ v on f2i, we have time-average total variation convergence to p, 
i.e., 




(b) In addition, if there exists no integer m > 1 such that [a, b] C [i, ^^j], (hence in partic- 
ular if a = 0), then we have convergence in total variation to p for every initial distribution 
¥ v on Oi ; i.e., 

lim sup \Vr,(rj(t) e A) - p{A)\ = 0. 



Proof. We prove this theorem by constructing a coupling. The two processes to be coupled 
have initial configurations r/ 1 and r/ 2 , with r] 1 ,^ 2 € fix. We denote by rj 1 ^), rj 2 (t) two 
independent copies of the process starting from rj 1 and rj 2 respectively. The corresponding 
independent additions at each time step are denoted by U} and U 2 , respectively. Let T\ = 
min{t : 77 1 (t) = 0} and T2 = min{t : r] 2 (t) = 0}. Suppose (without loss of generality) that 
T2 > T\. We define a shift-coupling ([H], Chapter 5) as follows: 



fj 1 ^) = Tj (t) for all t, 



v 2 (t) 




for t < T 2 , 
(Ta-Ti)) fort>T 2 . 



Defining T = min{t : rj 1 ^) = n 2 (t) = 0}, we also define the exact coupling 

^(t) = ^(t) for all t, 
( r] 2 (t) for t < T, 
} ri l {t) for t > T. 



fj 3 (t) 



Since the process is Markov, both couplings have the correct distribution. We write P = 
x P^ 2 . Since P(T 2 < 00) = P(Ti < 00) = 1, the shift-coupling is always successful, and 
(a) follows. 

To investigate whether r/ 1 (f) = rj 2 {t) = occurs infinitely often P-a.s., we define N = {n : 
(n — l)a < 1 A nb > 1}; this is the set of possible numbers of time steps between successive 
events ^(t) = 0. In words, an n € J\f is such that, starting from ry 1 = 0, it is possible that 
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in n — 1 steps we do not yet reach energy 1, but in n steps we do. To give an example, if 
a > 1/2, then TV = {2}. 

If the gcd of TV is 1 (this is in particular the case if a = 0), then the processes {t : 
ry 1 (t) = 0} and {t : r] 2 (t) = 0} are independent aperiodic renewal processes, and it follows 
that 77 1 (i) = 7] 2 (t) = happens infinitely often P-a.s. 

As we have seen, for a > 0, the gcd of TV need not be 1. In fact, we can see from the 
definition of TV that this is the case if (and only if) there is an integer m > 1 such that 
[a, b] C [— , —33;]. Then TV = {m}. For such values of a and b, the processes {t : ^(t) = 0} 
and {t : rj 2 (t) = 0} are periodic, so that we do not have a successful exact coupling. □ 

4.2 The stationary distribution of the (1, [0, b])-model 

We write p b for the stationary measure p ob of the (1, [0, 6])-model and F b for the distribution 
function of the amount of energy at stationarity, that is, 

F b (h) =p b (r] :0<r/<h). 

We prove the following explicit solution for F b (h). 

Theorem 4.2. (a) The distribution function of the energy in the (1, [0,b]) -model at station- 
arity is given by 



F\h) = < 
where mh = |~f ] — 1 and where 



forh< 0, 
F 6 (0) > forh = 0, 
F b (0)^ o t^(h-nbre^ for0<h<l, 

1 for h> 1, 



(17) 



^(0) = — — ^ 



follows from the identity F b (l) = 1. 
(b) For h € [0, 1] we have 



b 



lim F b (h) = h. 
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We remark that although in (a) we have a more or less explicit expression for F b (h), the 
convergence in (b) is not proved analytically, but rather probabilistically. 

Proof of Theorem part (a). Observe that the process for one site and a = is defined 
as 

rj(t + 1) = (rj(t) + Ut+i) lj 7 (t)+c/ t+1 <i- (18) 

We define F b = ¥(rj(t) < h), and derive an expression for F b +1 (h) in terms of F b (h). In 
the stationary situation, these two functions should be equal. We deduce from (fT8"j) that for 
< h< 1, 

F b +1 (h) = F( V (t) + Ut+i <h) + P(r?(t) + U t+1 > 1). (19) 
We compute for < h < b, 

P(r/(t) + U t+ i < h) = ¥( V (t) <h- U t+l ) 



h 1 



-F(r](t) <h-u)du 



£±F b (h-u)du, (20) 



and likewise for b < h < 1 we find 



F(r,(t) + U t+1 <h)= f \{F b {h - u)) du. (21) 
Jo o 

Finally, 

P««) + t/ t+ i>l) = I i(f*(l)-f*(l-«))«fo 

- jf (1 ~ f? t (1 "" )) *-^c)- w 



Putting (fTUj) . (|2"0"j) . (|2~Tj) and (|2"2~j) together leads to the conclusion that the stationary 
distribution F b (h) satisfies 



( j " *, + **(()) if&<*<l. 

Furthermore, since F b (h) is a distribution function, F b (h) = for h < and .F 6 (l) = 1. 
We can rewrite equation ([23]) as a differential delay equation. We take f b (h) = dF d ^ the 
density corresponding to F b for < h < 1; this density exists according to Theorem [ 
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We first consider the case < h < b, in which case mt = 0. We differentiate (f23j) twice, 
to get 

which leads to the conclusion that F b (h) = F b (0)e~b , consistent with (|17p for < h < b. 
Now we consider the case b < h < 1. We differentiate (|23p on both sides to get 

/ 6 (/ l ) = i(F 6 (/ l )-F fe (/ l -6)). (24) 

At this point, we can conclude that the solution is unique and could in principle be found 
using the method of steps. However, since we already have the candidate solution given in 
Theorem 14.21 we only need to check that it indeed satisfies equation (|23|) . We check that for 
the derivative f b of F b as defined in (fTTl) . for b < h < 1, 



f\h) 



F b (0) (<h - nbY^e^ + (h- KbY l -e^\ 

K = K ' ^ ' 

F b (0) ^, l._ ! 1 , L ink i 

F b (0) ^. l u l ,. K ^ 

+ 6 E(-ft) kiS h -^ e " ' 

K=0 



whereas 



and 



F b (h) F b (0) ^k. l. K 1 ; 

-r i = ^ i E(-^(^-^ 

K = 



K=0 

which leads to ([23]) as required. □ 

We remark that the probability density function f b (h) has an essential point of discontinuity 
at h = b. Figures [TJ and [2] show two examples of f b (h). 
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0.2 0.4 0.6 0.8 1 

Figure 1: f b (h) for b = ^. Note the discontinuity at h = 




0.2 0.4 0.6 0.8 1 



Figure 2: for b = jq. This figure illustrates that for small b, f b (h) tends to the uniform 

distribution. 

Proof of Theorem part (b). Without loss of generality, we start at time with zero 
energy. To avoid confusion and express the dependence of the models on the parameter b, 
we will write r) b (t) for the random state of the (1, [0, 6])-model at time t, with rj b = 0. 

In this proof, it is helpful to see that the process can be viewed as an alternating renewal 
process. To this end, we define the following random times, with i = 1,2,...: 

T b (0) := mm{n > 7^(0), v\ri) = 0}, 

with 7^(0) = 0. We also define 

T b {h) := max{n G (T b (0), T b +l (0)) : V b (s) < h,Vs G {T b (0) + 1, . . . ,n}}. 

We fix h. The process alternates between states where rj b (t) < h and r/ b (t) > h, such that 
renewal events occur at times Zf (0). All intervals between successive T b (0) are i.i.d., as they 
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are all intervals of the form T b (h) — Zf (0). Thus, the requirements of an alternating renewal 
process are met, and we can conclude that 

F b (h) = lim FA V b (t) <h) = ^MM>±1 (25) 

v ; *-»°o " w v ' - ' E(r|(o)) 

which is valid for all b < 1. 

To compute the expectations in (|25[) in the limit 6 — > 0, we use another process: let 
U\,U>2,... be independent uniform random variables on the interval [0,1] and write (as 
before) = U x + ■ ■ ■ + U n . We define 

N(s) = max{n G N : < s}. 

The process {N(s) : s > 0} is a renewal process and by the elementary renewal theorem, 

lim « =E([/l) 4 (26) 

Observe that T?(h) and iV(£) have the same distribution; this is just a rescaling. Likewise, 
r|(0) has the same distribution as -/V(i) + 1- Hence, 

Fb(h) _ + 1 (27] 

and by (|26p we find that 

lim — i — = 1. 

4 



Combining this with (j27[) we conclude that 

lim h 

6-+0 (E(N(l)) +1) -2b 



limF ^, _ limft .(«(iv(f)) + i)-f = ft . 



□ 



5 The (iV, [a, 6])-model with N > 2 and a > | 
5.1 Uniqueness of stationary distribution 

In the course of the process of Zhang's model, the energies of all sites can be randomly aug- 
mented through additions, and randomly redistributed among other sites through avalanches. 
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Thus at time t, every site contains some linear combination of all additions up to time t, 
and the energies at time 0. In a series of lemma's we derive very detailed properties of these 
combinations in the case a > 1/2. These properties are crucial to prove the following result. 

Theorem 5.1. The (N, [a,b]) model with a > \, has a unique stationary distribution fx = 
[i ab . For every initial distribution v on Ojy, F u converges exponentially fast in total variation 
to fj,. 

We have demonstrated for the case a > \ that after a finite (random) time, we only 
encounter regular configurations (Proposition 13.5] ). By Lemma 13.61 if i]{t — 1) is regular, 
then the knowledge of TZ(n(t — 1)) and Xt suffices to know the number of topplings of each 
site at time t. Thus, also the factors F^ in Proposition 13.71 are functions of lZ(n(t — 1)) and 
Xt only. Using this observation, we prove the following. 

Lemma 5.2. Let a > ^. Suppose at some (possibly random) time r we have a configuration 
£(r) with no anomalous sites. Then for all j = 1, . . . , N and for t > t, we can write 

t N 

= A 0^) U + E B mj(t)U(r) (28) 

9=T+l m=l 

in such a way that the coefficients in (J28]) satisfy 

A ej (t) = A ej (n(C(r)),X T+1 ,...,X t ) 

and 

B mj (t) = B mj (n{^T)),X T+1 ,. ..,X t ), 
and such that for every j and every t > r, 

t N 

e= T +i m=l 



Remark 5.3. Notice that, in the special case that r is a stopping time, Agj(t) is independent 
of the amounts added after time r, i.e., A$j(t) and {Uq, 9 > r + 1} are independent. We will 
make use of this observation in Section 15. 2[ 
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Proof. We use induction. We start at t = r, where we choose B m j(r) = T^-(C(, T ))j ^-m=j- We 
then have Ylm=i ^mjij) = 7Z(£,{ T ))j-, so that at t = r the statement in the lemma is true. 
We next show that if the statement in the lemma is true at time t > r, then it is also true 
at time t + l. 

At time t we have for every j = 1, . . . , N , 

t N 



&(t) = + E (r) 



=r+l 



m=l 



with 



^ ^-(t) + X) B ^(<) = ^(*))i. 

9=t+1 m=l 

where all Agj(t) and B m j(t) are determined by TZ(^(t)), X T+ i, . . . , X tl so that 7£(£(t)) is 
also determined by 7£(£(t)), X r+ i, . . . ,X t . We first discuss the case where we added to a full 
site, so that an avalanche is started. In that case, the knowledge of 7£(£(t)), X t+ i, . . . , Xt+i 
determines the sets C(t + 1), T(t + 1) and the factors from Proposition 13.71 (This last 
fact is due to the fact that a > 1/2.) We write, denoting X t +\ = x, 



E FiMt) + FxjUt+i 
ieT(m) 

" t N 

E * 

ieT(t+i) 



m=l 



+ FxjUt+i- 



Thus we can identify 



A ej (t + 1)= Y F v A ei(t), 
ieT(t+l) 

B m j(t + 1)= ^2 FijB m .i(t), 
ieT(m) 



(29) 
(30) 



and 



A t+1J (t + l) = F xj , 

so that indeed all Agj(t + 1) and B m j(t + 1) are functions of 7£(£(t)), X t+ i, . . . , 



X t+ i only. 
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Furthermore, 



t+i N 



r t N 



Fij YAgi^ + ^B^it) + F, 



xj 



0=1 m=l 



ieT(t+l) 18=1 m=l 



E F ij + F xj =n{r ]j {t + l)), 



ieT(t+l) 



(31) 



where we used that all sites that toppled must have been full, therefore had reduced value 
1. 

If no avalanche was started, then the only site that changed is the addition site x, and it 
must have been empty at time t. Therefore, we have for all 9 < t+1, Ag x (t+1) = Ag x (t) = 0, 
for all m, B mx (t + 1) = B mx {t) = and A t +i tX (t + 1) = 1, so that the above conclusion is 
the same. □ 

For every 9, we have Yl%eT(t+i) Aei(t) < 1, because the addition U$ gets redistributed by 
avalanches, but some part disappears through topplings of boundary sites. One might expect 
that, as an addition gets redistributed multiple times and many times some parts disappear 
at the boundary, the entire addition eventually disappears, and similarly for the energies 
£j(r). Indeed, we have the following results about the behaviour of Ag^t) for fixed 9, and 
about the behaviour of B m j(t) for fixed m. 

Lemma 5.4. For every 9, and for t > 9, 

1. maxi<j<Ar Agi(t) and maxi<j< n B m i(t) are both non-increasing int. 

2. For all 9, m and i, lim^oo Ag^t) = 0, and lim^oo B m i(t) = 0. 

Proof. We can assume that t > 9. The proofs for Agj(t) and for B m j(t) proceed along 
the same line, so we will only discuss Agj(t). We will show that for every j, Agj(t + 1) < 
maxj Agi(t), by considering one fixed j. If the energy of site j did not change in an avalanche 
at time t + 1, then 

A ej {t + 1) = A 6j (t) < max A ei (t). 

i 

If site j became empty in the avalanche, then 

Agj(t + 1) = < max A ei (t). 
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For the third possibility - the energy of site j changed to a nonzero value in an avalanche at 
time t + 1 - we use (|30|) , and estimate 

A ej (t + 1)= FijA 6i (t) <TnaxA ei (t) ^ 

ieT(t+i) ' i£T(t+l) 

By Proposition Opart (1) and (2), EieT(t+i) ^ - 1 ~ 2~ T3iv/2n ^ go that in thig third cag6) 
AgJt + l) < (1 - 2~ r37V/2l )max^(t) < max4(t). 

i i 

Thus, it follows that maxj^4gj(t + 1) can never be larger than maxj^4gj(t). This proves 
part (1). It also follows that when between t and t + r all sites have changed at least once, 
we are sure that maxj Agi(t + r) < (1 — 2~^ 3iV / 2 l ) maxj Agi(t). 

We next derive an upper bound for the time that one of the sites can remain unchanged. 
Suppose at some finite time t when r](t) is regular (see Proposition 13, 5|) . we try never to 
change some sites again. If all sites are full, then this is impossible: the next avalanche will 
change all sites. If there is an empty site x, then the next addition changes either the sites 
1, . . . ,x, or the sites x, . . . , N. In the first case, after the avalanche we have a new empty 
site x' < x. If we keep trying not to change the sites x, . . . , N, we have to keep making 
additions that move the empty site closer to the boundary (site 1). It will therefore reach 
the boundary in at most N — 1 time steps. Then we have no choice but to change all sites: 
we can either add to the empty site and obtain the full configuration, so that with the next 
addition all sites will change, or add to any other site, which immediately changes all sites. 
This argument shows that the largest possible number of time steps between changing all 
sites is N + 1. We therefore have 

max A ei (t) < (1 - 2-r 3Ar / 2 l)LifTrJ ; (32) 

i 

so that 

lim max A ei (t) < lim (1 - 2~r 3Ar / 2 l ) L^rJ = o. 

t^oo i t— >oo 

□ 

With the above results, we can now prove uniqueness of the stationary distribution. 

Proof of Theorem \5.1[ By compactness, there is at least one stationary measure [i. To 
prove the theorem, we will show that there is a coupling (r) 1 (i ) , ?7 2 (t))o° with probability law 
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P^i^) for two realizations of the (N, [a, b]) model with a > ^, such that for all e > 0, and 
for all starting configurations r/ 1 and r) 2 , for t — > oo we have 

P (r)W) (max - > = 0( e -°^)), (33) 

with a^v > 0. 

From (|33p . it follows that the Wasserstein distance ([3], Chapter 11.8) between any 
two measures /ii and /j.2 on fijv vanishes exponentially fast as t ^ 00. If we choose r/ 1 
distributed according to /i stationary, then it is clear that every other measure on fijv 
converges exponentially fast to fj,. In particular, it follows that /i is unique. 

As in the proof of Theorem l4.lt the two processes to be coupled have initial configurations 
rj 1 and 7] 2 , with r/ 1 ,^ 2 G Ojy. The independent additions at each time step are denoted by 
U} and U 2 , the addition sites X\ and X 2 . 

We define the coupling as follows: 



f) 1 (t) =r ? 1 (t) for all t 

J f] 2 {t) for f < T, 

{ A uhx i(fj 2 (t-l)) fort>T, 



where T = min{i > T' : IZ^ 1 ^)) = 1Z(rj 2 (t))}, and T' the first time that both 77 1 (i) and 
r] 2 (t) are regular. In Proposition 13.51 it was proven that T" < N(N — 1), uniformly in r]. In 
words, this coupling is such that from the first time on where the reductions of ?) 1 (i) and 
fj 2 (t) are the same, we make additions to both copies in the same manner, i.e., we add the 
same amounts at the same location to both copies. Then, by Lemma [3.61 in both copies the 
same avalanches will occur. We will then use Lemma 15.41 to show that, from time T on, the 
difference between r) 1 (t) and fj 2 (t) vanishes exponentially fast. 

First we show that P^i^2\(r > t) is exponentially decreasing in t. There are N + 
1 possible reduced regular configurations. Once fj (t) is regular, the addition site X^ +1 
uniquely determines the new reduced regular configuration 1Z{j] l (t + 1)). This new reduced 
configuration cannot be the same as 7?.(r/ 1 (i)). Thus, there are N equally likely possibilities 
for iKrj 1 ^ + 1)), and likewise for K(r?(t + 1)). 

If 7^.(7y 1 (i)) 7^ TZ(rj 2 (t)), then one of the possibilities for 7^(r/ 1 (i + 1)) is the same as 
7Z(r] 2 (t)), so that there are N — 1 possible reduced configurations that can be reached both 
from r? 1 (t) and rj 2 (t). The probability that 1Z{j] l (t + 1)) is one of these is ^r - , and the 
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probability that K(r] 2 (t + 1)) is the same is -j^. Therefore, T is geometrically distributed, 
with parameter pn = 

We now use Lemma 15.21 with r = T. For t > T, we have in this case that A l - e {t) = A 2 g(t) 
and Bj m (t) = B 2 m (t), because from time T on, in both processes the same avalanches occur. 
Also, for t > T, we have chosen U l (t) = U 2 (t). Therefore, for t > T, 

N 
m=l 

From (|32p in the proof of Lemma 15.41 we know that 
so that 

N 

^^(^^^(i-rWi)^, 

m=l 

so that for t > T, we arrive at 

max \fj)(t) - fj]{t)\ < 2N(1 - 2^ 3N ^)^r. 

We now split P( J) i J) 2)(maXj \fjj(t) — fjj(t)\ > e) into two terms, by conditioning on t < 2T 
and t > 2T respectively. Both terms decrease exponentially in t: the first term because 
the probability of t < 2T is exponentially decreasing in t, and the second term because for 
t > 2T, maxj \ry-(t) — fj 2 (t)\ itself is exponentially decreasing in t. □ 

A comparison of the two terms P(t < 2T) and maxj — yields that for large, 

the second term dominates. We find that qln depends, for large N, on N as = — \ ln(l — 
2 _ l" 3Ar / 2 l )tv+t. We see that as N increases, our bound on the speed of convergence decreases 
exponentially fast to zero. Furthermore, we remark that from the above proof it also follows 
that all stationary temporal correlations decay exponentially fast, i.e., there exist /3at > 
such that for all square integrable functions /(??), with f f 2 dv = 1, J fdu = 0, 

5.2 Emergence of quasi-units in the infinite volume limit 

In Proposition 13.21 we already noticed a close similarity between the stationary distribution 
of Zhang's model with a > 1/2, and the abelian sandpile model. We found that the stationary 
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distribution of the reduced Zhang's model, in which we label full sites as 1 and empty sites 
as 0, is equal to that of the abelian sandpile model (Proposition 13. 5|) , 

In this section, we find that in the limit N — ► oo, the similarity is even stronger. We find 
emergence of Zhang's quasi-units in the following sense: as N — ► oo, all one-site marginals 
of the stationary distribution concentrate on a single, nonrandom value. We believe that 
the same is true for a < 1/2 also (see Section I6T31 for a related result), but our proof is not 
applicable in this case, since it heavily depends on Proposition 13.51 To state and prove our 
result, we introduce the notation /zjy for the stationary distribution for the model on N sites, 
with expectation and variance E^ and Var^, respectively. 

Theorem 5.5. In the (N, [a,b]) model with a > |, for the unique stationary measure /i/y 
we have 

lim fi N = <W) (34) 

where (W;/) denotes the Dirac measure concentrating on the (infinite-volume) constant con- 
figuration rji = E(C7) for all i £ N, and where the limit is in the sense of weak convergence 
of probability measures. 

We will prove this theorem by showing that for r\ distributed according to /xjv, in the 
limit N — > oo, for every sequence 1 < jjy < N, 

1. lim N ^ 00 E N ri jN =EU, 

2. lim^oo Var iV (r? ijv ) = 0. 

The proof of the first item is not difficult. However, the proof of the second part is compli- 
cated, and is split up into several lemma's. 

Proof of Theorem \5. 51 part (1). We choose as initial configuration r] = 0, the configuration 
with all iV sites empty, so that according to Lemma 15.21 we can write 

ViJt) = Y,A ejN (t)U d . (35) 

6=1 

Denoting expectation for this process as E^, we find, using Remark 15.31 that 

E$ VjN (t)=W E$K( V (t)) jN . 
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First, we take the limit t — > oo. By Theorem 15.11 E,Qt]j N (t) converges to E N r]j N . From 
Proposition 13.51 it likewise follows that rimt_ >00 EQlZ(rj(t))j = Inserting these and 

subsequently taking the limit N — > oo proves the first part. □ 

For the proof of the second, more complicated part, we need a number of lemma's. First, 
we rewrite Var N (rjj N ) in the following manner. 

Lemma 5.6. 

Var N (rj jN ) = Varft/) lim Eq 

t— >oo 



+ (EUf 



N 



(N + l) 



Proof. We start from expression (|35p . and use that the corresponding variance Var^ con- 
verges to the stationary Var as t — > oo by Theorem 15.11 We rewrite, for fixed N and 
3N = j, 



Vatf^Ct)) = [( Vj (t)) 2 ] - [E^-(t)] 2 = 



N 



E, 



E(U 2 )E% 



Y^{A ej {t)f u e + E A*j(t)UoA*>i{t)Uo, 
e=i 8^8' 



E(^)) 5 



+ {EU) 2 E% 



8=1 

- (EC/) 2 



(E(C/ 2 ) - (EU) 2 ) E$ 



E(^)) 5 



(E?7) z 



0=1 



V 



Var(C7) Eg 



+ (EC7) 2 Var^(^(7 ? (t)) J )- 



Where in the third equality we used the independence of the ^.-coefficients of the added 
amounts Uq. We now insert j = jjy, take the limit t — > oo, and insert lim^oo Var^ (TZ(r/(t))j 



Var N (1l(7]) 



— N 
3NJ — (N+l) 2 ■ 



JJN) 
□ 
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Arrived at this point, in order to prove Theorem 15.51 it suffices to show that 



0. (36) 



hm limE^ J2(Ag jN (t)f 

iV— >oo t^oo z — » 

Le=i 

The next lemma's are needed to obtain an estimate for this expectation. We will adopt 
the strategy of showing that the factors Agj(t) are typically small, so that the energy of a 
typical site consists of many tiny fractions of additions. To make this precise, we start with 
considering one fixed > N(N — 1), a time t > 0, and we fix e > 0. 

Definition 5.7. We say that the event Gt(a) occurs, if max j Agj(t) < a. We say that the 
event i?t(e) occurs, if m&xj Agj(t) > e, and if in addition there is a lattice interval of size at 
most M = [=] + 1, containing Xg, such that for all sites j outside this interval, Agj(t) < e. 
(We call the mentioned interval the 6-heavy interval.) 

Note that since we have Agj(t) < 1 for every 0, the number of sites where A$j(t) > e, 
cannot exceed ["-]. In Lemma 15.41 we proved that maxj Agj(t) is nonincreasing in t, for 
t > 0. Therefore, also Gt(ct) is increasing in t. This is not true for Ht(e), because after an 
avalanche, the sites where Agj(t) > e might not form an appropriate interval around Xg. 

In view of what we want to prove, the events Gt(e) and H t (e) are good events, because 
they imply that (if we think of N as being much larger than M) Agj(t) < e 'with large 
probability'. In the case that Gt(e) occurs, Agi{t) < e for all i, and in the case that Ht(e) 
occurs, there can be sites that contain a large Agi(t), but these sites are in the 0-heavy 
interval containing Xg. This latter random variable is uniformly distributed on {1, . . . ,N}, 
so that there is a large probability that a particular j does not happen to be among them. 
If we only know that Gt{p) occurs for some a > e, then we cannot draw such a conclusion. 
However, we will see that this is rarely the case. 

Lemma 5.8. For every N , for every > N(N — 1), for every e > 0, for every K and j, 
1. there exists a constant c = c{e), such that for < t < + K 



>»(Ag j (t)>e)<^; 



cK 



2. for every N large enough, there exist constants w = w(e) and < 7 = j(N,e) < 1, 
such that for t > 

P^(^-(t)> e )<(l- 7 )*- e - 3 -. 
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In the proof of Lemma 15.81 we need the following technical lemma. 

Lemma 5.9. Consider a collection of real numbers yi > 0, indexed by N, with < 1 

and such that for some x £ N, max^. < a. Then, for j > x + |~^] , we /iawe 



i_:E+2 1 / 1 \ 

E ^+ 2 ^ /(a):= V 1_ ^Tj" 



Proof. We write 



-x+2 \-\ j-x+2 

1 1 j-aH-2 

E^-+ 2 + ^liTT E W-i 
i=i 2 U *=r=i+i 



(37) 



Note that index x is in the second sum. For i = 1, . . . , |"~] , write 

Uj-i+2 = a - Zi, 

with < Zi < a. Then, since a\^~\ > 1 and ^ j/j < 1, we have I^"^ 2 , 1 Vj-i+2 < X][=i ^ 
so that 



E^%-+ 2 + i[T^ E %-+ 2 ^E^( a -^) + ^iiTTE ; 
t=i - i=fi]+i »=i - »=i 



r 1 ! r-i r-i 

where in the last step we used Z{ > 0. Thus 

E ^E^=( v 1 -^tJ«- 

1 = 1 4=1 v ' 



□ 



35 



Proof of Lemma [578{. part (1). We first discuss the case t = 9. We show that Hg(e) occurs, 
for arbitrary e. 

At t = 8, the addition is made at Xg. From Proposition 13.71 part (3), it follows, for every 
e, that if after an avalanche there are sites j with Agj(9) > e (we will call such sites '0-heavy' 
sites), then these sites form a set of adjacent sites including Xg, except for a possible empty 
site among them. Since we have Ylj=i Agj(9) < 1, there can be at most 6>-heavy sites. If 
the addition was made to an empty site, then Agx g (9) = 1. Thus, the #-heavy interval has 
length at most 1 + |~~] , and we conclude that Hg(e) occurs. To estimate the probability that 
Agj{9) > e, or in other words, the probability that a given site j is in the #-heavy interval, we 
use that Xg is uniformly distributed on {1, . . . , N}. Site j can be in the 9— heavy interval if 
the distance between Xg and j is at most 1+ so that (A 9j (0) > e) < 2±±^ =: 2^2. 

We next discuss 8 < t < 9 + K . We introduce the following constants. We choose a 
number w such that f w (l) < e, with / as in Lemma 15.91 and where f w denotes / composed 
with itself w times. Note that this is possible because lim^oo f k (l) = 0. We choose a 
combination of eo and d such that e w < e, with eu+i defined as + ^^-(/^(l) — Finally, 
we define M k = [i] + 1 + k(l + d). 

For fixed 9 and a time t > 9, we define three types of avalanches at time t: 'good', 
'neutral' and 'bad'. 

Definition 5.10. For a fixed 9, the avalanche at time t is 

• a good avalanche if the following conditions are satisfied: 

1. X t and Xg are on the same side of the empty site (if present) at t — 1, 

2. Xg is at distance at least M w from the boundary, 

3. Xt is at distance at least w from the boundary, and from the empty site (if present) 
at t — 1, 

4- Xt is at distance at least M w + [-] from Xg, 

5. Xg is at distance at least M w from the empty site (if present) at t — 1, 

• a neutral avalanche if condition (5) is satisfied, but (1) is not, 

• a bad avalanche in all other cases. 

Having defined the three kinds of avalanches, we now claim the following: 



36 



• If i?t_i(e) occurs, then after a neutral avalanche, Ht(e) occurs. 

• If Ht-i(e) occurs, then after a good avalanche, Gt(e) occurs. 

The first claim (about the neutral avalanche) holds because if condition (5) is satisfied, 
but (1) is not, then not only is Xg not in the range of the avalanche, but the distance of Xg 
to the empty site is large enough to guarantee that the entire 0-heavy interval is not in the 
range. Thus, the #-heavy interval does not topple in the avalanche. It then automatically 
follows that Ht(e) occurs. 

To show the second claim (about the good avalanche), more work is required. We break 
the avalanche up into waves. Using a similar notation as in the proof of Proposition 13.71 we 
will denote by Agj(k) the fraction of Ug at site j after wave k. We also define another event: 
we say that H k (M k , d k , e k ) occurs, if max j Agj(k) < a(k), and all sites where AgAk) > e 
are in an interval of length at most M k containing Xg (we will call this the (£;-#)-heavy 
interval), with the exception of site Xt when it is unstable, in which case we require that 
Ag Xt (k) < 2e fc . 

We define a 'good' wave, as a wave in which all sites of the (k-6)-heavy interval topple, 
and the starting site X t is at a distance of at least ^- from the (fc-0)-heavy interval. It might 
become clear now that Definition 15.101 has been designed precisely so that a good avalanche 
is an avalanche that starts with at least w good waves. We will now show by induction in 
the number of waves that after an avalanche that starts with w good waves, Gf(e) occurs. 

For k = 0, we choose do = 1, so that at k = 0, Hq{Mq, do, e~o) occurs. We will choose 
dfc+i = /(dfc) and e k +i = £fc + -^r+ri&k ~~ ^k)i so that once H W (M W , 

occurred, we 

are sure that after the avalanche Gt(e) occurs, because both a w and e~ w are smaller than 
e. Now all we need to show is that, if Hk(Mk, df~, ejfc) occurs, then after a good wave 
H k+ i(M k+1 ,a k+1 ,e k+1 ) occurs. 

From ©, we see that, for all j > X t that topple (and do not become empty), 

Ag 3 {k + 1) = -Ae ij+ i(k) + \A dj {k) + -Ae^k) + ■■■+ 2J } Xt+2 A e>Xt {k), (38) 
and similarly for j < X t . For j = X t , we have 

A ej (k+1) = ^A e j +1 (k)+^A j>0 (k))l Ag . +imo + (39) 

First we use that in a good wave, all sites in the 0-heavy interval topple, and Xt is 
not in this interval. We denote by m the leftmost site of the #-heavy interval, so that the 
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rightmost site is m + M(k). Suppose, without loss of generality, that X% < m. We substitute 
Agj(k) < ce k for all j in the 0-heavy interval, Ax t g(k) < 2e k , and Agj(k) < e k otherwise into 
([38]) . to derive for all j that topple: 

j<m-l,j^x A 9j (k + l)<e k , 

j = m-l,...,m + M(k) A ej (k + 1) < a k , (40) 
j = m + M(k) + d' A ej (k + 1) <e k + ^(a k -e k ), d' = l,2,... 

Additionally, by ([39]) . we have Ag t xt(k + 1) < 2e k . The factor 2 is only there as long as 
site X t is unstable. From (f4T)l) we have that a k+ \ < a k , but moreover, in a good wave, 
the variables Agj(k) satisfy the conditions of Lemma Ell so that in fact a k+ \ < f(a k ). If 
we insert our choice of d for d' , then we can see from ()40|) that indeed after the good wave 
H k+ i(M k+1 ,a k+ i,e k+1 ) occurs. 

Now we are ready to evaluate {Agj{t) > e), for t £ {6 + 1, . . . , + K}. As is clear by 
now, there are three possibilities: Gt(e) occurs, so that maxj Agj(t) < e, or H t (e) occurs, in 
which case Agj(t) can be larger than e if j is in the 0-heavy interval. We derived in the case 
t = 9 that the probability for this is bounded above by ^. Finally, it is possible that neither 
occurs, in which case we do not have an estimate for the probability that A$j(t) > e. But 
for this last case, we must have had at least one bad avalanche between 6 + 1 and t. We 
will now show that the probability of this event is bounded above by -jf*, where c\ depends 
only on e. 

As stated in Definition 15.101 a bad avalanche can occur at time t if at least one of the 
conditions (2) through (5) is not satisfied. Thus, we can bound the total probability of a 
bad avalanche at time t, by summing the probabilities that the various conditions are not 
satisfied. We discuss the conditions one by one. 

• The probability that condition (2) is not satisfied, is bounded above by since Xq 
is distributed uniformly on {1, . . . , iV}. 

• The probability that condition (3) is not satisfied, is bounded above by since X t 
is distributed uniformly on {1, . . . , iV}, and independent of the position of the empty 
site at t — 1, if present. 

• The probability that condition (4) is not satisfied, is bounded above by 2 ^ Iw ^i^ ; 
since Xt and Xq are independent. 
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• The probability that condition (5) is not satisfied is bounded above by ^4f-, since the 
position of the empty site at t — 1 is uniform on {1, . . . , N}. 

Thus, the total probability of a bad avalanche at time t is bounded by + t^t + 2 (^"^M) + 
= ft, so that the probability of at least one bad avalanche between 0+1 and t is bounded 
by We conclude that for t 6 {9 + 1, . . . , 9 + K}, F$(A ej (t) > e) < ^ft±^ < f or 
some c > 0. □ 

Proof of Lemma \5.8\ part (2). From Lemma EU it follows that if Gt(a) occurs, and after s 
time steps all 0-heavy sites have toppled at least once, in avalanches that all start at least 
a distance \^] from all current 0-heavy sites, then Gt+ S (f(a)) occurs. We will exploit this 
fact as follows. 

Suppose that Gt(a) occurs and that in addition, at time t the distribution of the empty 
set - if present - is uniform on {1, . . . , N}. We claim that this implies that Gt+2{f {ot)) occurs 
with a probability that is bounded below, uniformly in N and e. To see this, observe that if 
there is no empty site, then all N sites topple in one time step. If there is an empty site, this 
(meaning all sites topple) also happens in two steps if in the first step, we add to one side 
of the empty site, and in the second step to the other side. Denote by e\ the position of the 
empty site before the first addition (at Xt+i), and e2 before the second addition (at X t +2)- 
If Xt+i < ei, then &i < e\. Therefore, all sites topple if Xt+i < e\ and Xt+2 > e±. With 
the distribution of e\ uniform on {1, . . . , iV}, the probability that this happens is bounded 
below by some constant 7' independent of iV and e. However we have the extra demand 
that both additions should start at least a distance \— ] < [-] from all current 0-heavy sites, 
of which there are at most [=]. Thus, both additions should avoid at most [^] 2 sites. The 
probability that this happens is therefore less than some 7' > 0, but it is easy to see that 
the difference decreases with N. We can then conclude that there is an N' large enough so 
that the probability that this happens is at least 7 > for all N > N', with < 7 < 1 
independent of N and e. 

In view of this, the probability that Gg + 2(f(l)) occurs, for N large enough, is at least 7. 
We wish to iterate this argument w times. However, the lower bound 7 is only valid when 
the distribution of the empty site, if present, is uniform on {1, ... , iV}. We do not have this 
for r](9 + 2) since we have information about what happened in the time interval (9, 9 + 2). 
However, after one more addition, the position of the empty site in rj{9 + 3), if present, is 
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again uniform on {1, . . . , JV}. Since f w (l) < e, iterating this argument gives 



>$(maxA ej (t)>e) < (1 - 7) 
3 



t-e-3w 



□ 



Proof of Theorem \5.5\, part (2). By Lemma I5U1 it suffices to prove ([36]) . We estimate, using 
that Ee A ej N (t) < 1, 



E 



N 



E^w*)) ! 



< E 



JV 



max ^e,- (t) 
t-K<e<t JN 



+ E, 



— 0=1 



< E, 



N 



max A ej (t) 

K6<t 



max 

N(N+l)<e<t-K 



max ^6»7jv(0 



We then for the first two terms estimate, using that maxg Agj N (t) < 1, 

E^ [maxA ejN {t)] < e + F$ (maxA 6jN (t) > e) < e + E< ( A o jN (t) > e). 

We finally use Lemma 15.81 and choose IT = Kn increasing with JV. For 6 € [f,t — if at], we 
straightforwardly obtain Y^e=t-K N ^0 (-^VjjvC*) > e ) = 0(.~w~)> uniformly in t, as JV — > 00. 
For 6 < t — if/v we calculate 



E P S r (^iir(*)>e)< E (l-7)'^=0((l- 7 )^), 
6<t-K N t-9>K N 



JV -> 00, 



so that 



E, 



E^W')) S 



=1 



< 2e+ 0(^)+0((l- 7 )^)+E^ 



max ^^(t) 



JV 00. 



In the limit t — > 00, by Lemma 15.41 part (2), the last term vanishes. We now choose 
iiQv = JV 1 / 3 , to obtain 



A' 



limsup lim E 

AT-KX3 t^oo 



E( A ^(*)) 2 



Since e > is arbitrary, we finally conclude that 



lim lim E^ 

N^oo t^oo 



E^(*)) S 



< 2e, 



0. 



□ 
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6 The (N, [0, l])-model 



6.1 Uniqueness of the stationary distribution 

Theorem 6.1. The (A, [0, 1]) model has a unique stationary distribution vn- For every 
initial distribution v on Qjy, P u converges in total variation to vjy. 

Proof. We prove this theorem again by constructing a successful coupling. For clarity, we 
first treat the case A = 2, and then generalize to A > 2. The coupling is best described in 
words. 

Using the same notation as in previous couplings, we call two independent copies of the 
process rfit) and rj 2 (t), and call the coupled processes fj 1 ^) and fj 2 (t). Initially, we choose 
f} 1 (t) = r/ 1 (t), and r) 2 (t) = r] 2 (t). It is easy, but tedious, to show that ri\(t) = i] 2 (t) = 0, 
while lZ{n\{t)) = 1Z(i] 2 {t)) = 1, occurs infinitely often. 

At the first such time T\ that this occurs, we choose the next addition as follows. Call 
A(t) = r)l(t) - r)l{t). We choose X% 1+1 = X* i+1 , and [>| 1+1 = (U^_ 1+1 + A(Ti))mod 1. 
Observe that the distribution of U^+i 1S uniform on [0, 1]. 

This addition is such that with positive probability the full sites are chosen for the 
addition, and the difference A(Ti) is canceled. More precisely, this occurs if X^ 1+1 = 2, 
which has probability 1/2, and (C/j. +1 +A(Ti))mod 1 = f/j. +1 +A(Ti), which has probability 
at least 1/2, since r]\(Ti) and i] 2 (Ti) are both full, therefore A(Ti) < 1/2. If this occurs, 
then we achieve success, i.e., 7? 1 (Ti + 1) = i? 2 (?i + 1), and from that time on we can let the 
two coupled processes evolve together. 

If 57 1 (Ti + 1) 7^ ffiTi + 1), then we evolve the two coupled processes independently, and 
repeat the above procedure at the next instant that f)l(t) = fj 2 (t) = 0. Since at every such 
instant, the probability of success is positive, we only need a finite number of attempts. 
Therefore, the above constructed coupling is successful, and this proves the claim for N = 2. 

We now describe the coupling in the case N > 2. We will again evolve two processes 
independently, until a time where n\(t) = r] 2 (t) = 0, while all other sites are full. At this 
time we will attempt to cancel the differences on the other — 1 sites one by one. We 
define Aj(t) = r]j(t) — r] 2 (t), and as before we would be successful if we could cancel all 
these differences. However, now that N > 2, we do not want an avalanche to occur during 
this equalizing procedure, because we need r]\(t) = r] 2 (t) = during the entire procedure. 
Therefore, we specify T\ further: T\ is the first time where not only r]\{t) = rj 2 (t) = and 
all other sites are full, but also r/j(i) < 1 — e and rj 2 (t) < 1 — e, for all j = 2, . . . , A, with 
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e = 2 n+i ■ At such a time, a positive amount can be added to each site without starting an 
avalanche. We will first show that this occurs infinitely often, which also settles the case 
N = 2. 

By Proposition 13.21 after a finite time r/ 1 (t) and rj 2 (t) contain at most one non-full site. 
It now suffices to show that for any E Qjy with at most one non-full site, with positive 
probability the event that £i(t + 4) = 0, while £j(t + 4) < 1 — 2 -^ TT for every 2 < j < N, 
occurs. 

One explicit possibility is as follows. The first addition should cause an avalanche. This 
will ensure that + contains one empty site. This occurs if the addition site is a full site, 
and the addition is at least 1/2. The probability of this is at least |(1 — jj). The second 
addition should change the empty site into full. For this to occur, the addition should be 
at least 1/2, and the empty site should be chosen. This has probability to. The third 
addition should be at least 1/2 to site 1, so that an avalanche is started that will result in 
£jv(t + 3) = 0. This has again probability to. Finally, the last addition should be an amount 
in [^, §], to site N — 1. Then by ([!]), every site but site N will topple once, and after this 
avalanche, site 1 will be empty, while every other site contains at most 1 — 2 n+i • This last 
addition has probability to. 

Now we show that at time T% defined as above, there is a positive probability of success. 
To choose all full sites one by one, we require, first, for all j = 2, . . . , N that Xj, 1+ j_i = j- 
This has probability (jy)^ -1 . Second, we need (C4 1+J _ 1 +A i (Ti))mod 1 = U^+^+Ajpi) 
for all j = 2, . . . , N. This event is independent of the previous event and has probability at 
least {\) N ~ l - If this second condition is met, then third, we need to avoid avalanches, so for 
all j = 2, . . . , N, 7?)(Ti +j-l) + C4 1+i _i = Vj{Ti + j - 1) + U^+j-i < L lt is not hard 
to see that this has positive conditional probability, given the previous events. We conclude 
that the probability of success at time T\ + N — 1 is positive, so that we only need a finite 
number of such attempts. Therefore, the coupling is successful, and we are done. □ 

6.2 Simulations 

We performed Monte Carlo simulations of the (N, [0, l])-model, for various values of N. 
Figure [3] shows histograms of the energies that a site assumes during all the iterations. We 
started from the empty configuration, but omitted the first 10% of the observations to avoid 
recording transient behavior. Further increasing this percentage, or the number of iterations, 
had no visible influence on the results. 
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Figure 3: Simulation results for the (iV, [0, l])-model. The histograms represent observed 
energies during 100,000 (a,b) and 200,000 iterations (c-f). The system size was 3 sites (a,b), 
30 sites (c,d) and 100 sites (e,f). (a),(c) and (e) are boundary sites, (b), (d) and (f) are 
central sites. 

The presented results show that, as the number of sites of the model increases, the energy 
becomes more and more concentrated around a value close to 0.7. In the next section, we 
present an argument for this value to be y/l/2. We further observe that it seems to make a 
difference where the site is located: at the boundary the variance seems to be larger than in 
the middle. 

6.3 The expected stationary energy per site as iV — > oo 

From the simulations it appears that for large values of N, the energy per site concentrates 
at a value close to 0.7, for every site. Below we argue, under some assumptions that are 
consistent with our simulations, that this value should be y/l/2. 

First, we assume that every site has the same expected stationary energy. Moreover, 
we assume that pairs of sites are asymptotically independent, i.e., r\ x becomes independent 
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of r\ y as \x — y\ — > oo. (If the stationary measure is indeed such that the energy of every 
site is a.s. equal to a constant, then this second assumption is clearly true.) With E VN 
denoting expectation with respect to the stationary distribution ujv, we say that (vn)n is 
asymptotically independent if for any 1 < xn-,Un < N with \xn — yjy\ — » oo, and for any 
A, i? subsets of R with positive Lebesgue measure, we have 



lirn^ [E VN (1 Vxn eB \rj yN €A)- E VN (l„ Xjv e ^J = 0. (41) 

Theorem 6.2. Suppose that in the (N, [0, 1]) model, for any sequence jn € {1, . . . , N}, 

lim E VN (r, jN ) = p, (42) 

for some constant p. Suppose in addition that (vn)n is asymptotically independent. Then 
we have p = y^. 

Proof. The proof is based on a conservation argument. If we pick a configuration according 
to vm and we make an addition U, we denote the random amount that leaves the system by 
E ut,N- By stationarity, the expectation of U must be the same as the expectation of E outj N- 

The amount of energy that leaves the system in case of an avalanche, depends on whether 
or not one of the sites is empty (or behaves as empty) . Remember (Proposition 13. 2D that 
when we pick a configuration according to the stationary distribution, then there can be at 
most one empty or anomalous site. If there is one empty site, then the avalanche reaches 
one boundary. If there are only full sites, then the avalanche reaches both boundaries, and 
in case of one anomalous site, both can happen. 

However, configurations with no empty site have vanishing probability as N — > oo: we 
claim that the stationary probability for a configuration to have no empty site, is bounded 
above by pn, with limTv^ooPTV = 0. To see this, we divide the support of the stationary 
distribution into two sets: £, the set of configurations with one empty site, and M, the set of 
configurations with no empty site. The only way to reach TV from £, is to make an addition 
precisely at the empty site. As X is uniformly distributed on {1, ... , N}, this has probability 
i, irrespective of the details of the configuration. The only way to reach £ from M, is to 
cause an avalanche; this certainly happens if an addition of at least 1/2 is made to a full 
site. Again, since X is uniformly distributed on {1, . . . , N}, and since there is at most one 
non-full site, this has probability at least | N ^ . 

Now let X be the (random) addition site at a given time, and denote by A x the event 
that X = x and that this addition causes the start of an avalanche. Since E ou t.N = when 
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no avalanche is started, we can write 



N 



^v N (Eout,N) = ^ E VN (Eout,N\A x )P VN (Ac). 
x=l 

We calculate P VN (A x) as follows, writing U for the value of the addition: 
IPVijv {Ax) 



(43) 



^P^fe + U > 1) = ^P UJV (C/ > 1 - Vx) 



1 



1 



— J W VN (U > 1 - 7] x )dv N (r]) = — I ri x dv N (r]) 
—E VN (rj x ). 



(44) 



Let Ln = [log AT]. Even if the avalanche reaches both boundary sites, the amount of energy 
that leaves the system can never exceed 2, which implies that 



N 



N-2L 



N 



It follows from (14*3 



^2^v N (Eout,N\A x ) — ^2 ^v N {Eout,N\A x 
x=l x=2Lpf 

and that 

N-2L N 



< 8Ljv. 



(45) 



TV -2^ 

^(^.jv) = -r= E„ JV (^ )iV |A^)E„ iV (7 7:c ) + 0(L J v/iV). 



(46) 



x=2L A 



If the avalanche, started at site x, reaches the boundary at site 1, then the amount 
of energy that leaves the system is given by \r\\ + \r\2 + • • • + (^) x (Vx + U). For all x € 
{2Ljv, . . . , N — 2L/v}, this can be written as 

\m + \m + --- + (\) Ln vl n + (\) Ln+1 vl n+ i + ■■■ + {\f{rix + U), 

where for the last part of this expression, we have the bound 

{\) LN+X r}L N+ x + ■■■ + (\f(rix + U)< (\) L *. 

Since the occurrence of A x depends only on rj x (and on X and U), for 2L^ < x < N—2Ln, by 
asymptotic independence there is an c^TVi with lirnjv_ >00 ccn = 0, such that for all 1 < i < Ln 
and 2Ljv < x < N — 2Ln, we have 



E VN (r]i\A x ) -E VN (r]i)\ < a N , 
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so that 

(\m + ■■■ + (\) LN r]L N \A x ) - E VN (± m + ■■■ + (\) Ln Vl, 



which is bounded above by a at. By symmetry, we have a similar result in the case that the 
other boundary is reached. In case both boundaries are reached, we simply use that the 
amount of energy that leaves the system is bounded above by 2. 
In view of this, we continue the bound in (|46p as follows 

N-2L N /ii i i \ 

E VN (E out:N ) = - Y, (^v n (^Vi + 1 V2 + --- + (^ Ln Vl n \A x )+0((-) Ln ))^v n (.V x ) + 

x=2L N ^ ' 

+0(L N /N) 

1 N ~ 2Ln 1 1 1 
= N 12 E v N (^ + jV2 + --- + (^) LN 7]L N )E VN (r ]x ) + 

x=2Lpj 

+0(^) + 0({ l -) L ") + 0(a N ) + 0(p N ), 



as N — > oo. Letting N — > oo and inserting (|42p now gives 

lim (Eout.N*) = P 2 - 

As the expectation of U is i, we conclude that p = \/h- D 
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